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Abstract 

In the teleparallel equivalent of general relativity the energy density of asymptot- 
icaly flat gravitational fields can be naturally defined as a scalar density restricted 
to a three dimensional spacelike hypersurface E. Integration over the whole S yields 
the standard ADM energy. Here we obtain the formal expression of the localized 
energy for a Kerr black hole. The expression of the energy inside a surface of con- 
stant radius can be explicitly calculated in the limit of small a, the specific angular 
momentum. Such expression turns out to be exactly the same as the one obtained 
by means of the method proposed recently by Brown and York. We also calculate 
the energy contained within the outer horizon of the black hole, for any value of 
a. The result is practically indistinguishable from E = 2Mj r , where Mj r is the 
irreducible mass of the black hole. 
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I. Introduction 

Although it is widely believed that Einstein's equations describe the dynamics of the 
gravitational field, it has not been possible so far to arrive at a definite expression for the 
gravitational energy in the context of Einstein's general relativity. Attempts based on 
the Hilbert-Einstein action integral fail to yield an expression for the gravitational energy 
density^ 0. The total gravitational energy is normally obtained from surface terms in 
the action or in the Hamiltonian[]3], |J, or from pseudotensor methods which make use of 
coordinate dependent expressions. 

Recently an expression for quasi- local energy has been proposed by Brown and York||. 
Such expression is derived directly from the action functional A c \. The latter is identified 
as Hamilton's principal function and, in similarity with the classical Hamilton- Jacobi 
equation, which expresses the energy of a classical solution as minus the time rate of 
the change of the action, the quasilocal gravitational energy is identified as minus the 
proper time rate of change of the Hilbert-Einstein action (with surface terms included). 
Expressions for the quasilocal energy have been obtained for the Schwarzschild solution || 
and for the Kerr solution 0. 

Einstein's equations can also be obtained from the teleparallel equivalent of general 
relativity (TEGR). The Lagrangian formulation of the TEGR is established by means of 
the tetrad field e a M and the spin affine connection u^ a b, which are taken to be completely 
independent field variables, even at the level of field equations. This formulation has 
been investigated in the past in the context of Poincare gauge theories]?], ||. However, 
as we will explain ahead, this is not an alternative theory of gravity. This is just an 
alternative formulation of general relativity, in which the curvature tensor constructed 
out of u^ a b vanishes, but the torsion tensor is non-vanishing. The physical content of the 
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theory is dictated by Einsten's equations. In this alternative geometrical formulation the 
gravitational energy density can be naturally defined. 

The expression for the gravitational energy density arises in the framework of the 
Hamiltonian formulation of the TEGR|§. It has been demonstrated that under a suitable 
gauge fixing of ui^, already at the Lagrangian level, the Hamiltonian formulation of the 
TEGR is well defined. The resulting constraints are first class constraints ||. The Hamil- 
tonian formulation turns out to be very much similar to the usual ADM formulation [|3|]. 
However there are crucial differences. The integral form of the Hamiltonian constraint 
equation C = in the TEGR can be written in the form C = H — Eadm — 0, when 
we restrict considerations to asymptotically flat spacetimes||10||. The quantity e(x) which 
appears in the expression of C and which under integration yields Eadm is recognized as 
the gravitational energy density. 

We have calculated the energy inside a sphere of radius r a in a Schwarzschild spacetime 
by means of e{x) ||10|| . The expression turns out to be exactly the same as the one obtained 
by means of the procedure of ref.|| (expression (6.14) of ||). In this paper we consider the 
Kerr black-hole. We obtain the formal expression for the energy contained in any space 
volume in terms of non-trivial integrals. In the limit of slow rotation (small specific angular 
momentum) the energy contained within a surface of constant radius r a can be calculated. 
Again the result obtained here is exactly the same as that obtained by Martinez who 
adopted Brown and York's procedure. The advantage of our procedure rests on the fact 
that the localized energy associated with a Kerr spacetime can be calculated in the general 
case, without recourse to particular limits, at least by means of numerical integration, 
whereas in Brown and York's procedure one has to calculate the subtraction term e° and 
for this purpose it is necessary to embed an arbitrary two dimensional boundary surface 
of the Kerr space S in the appropriate reference space (E 3 , say), which is not always 
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possible 0. 

We have also calculated the energy contained within the outer horizon of the black 
hole. Such a quantity has been obtained by Martinez |J in the limit of small a, and reads 
E = 2Mi r (plus corrections of order 0(xhr) ), where Mj r is the irreducible mass of the 



black hole. The concept of irreducible mass was introduced by Christodoulou|[Tl|| . He 
showed that the mass of a rotating black hole cannot be decreased to values below M ir by 
means of Penrose's process of extraction of energy. One would thus consider E = 2M ir 
to be the energy that cannot escape from the black hole. Here we obtain the expression 
of the energy contained within the horizon for any value of a. The result is striking. The 
numerical values of this expression are practically coincident with 2Mj r in the whole range 
< a < m, although the expression is is algebrically different from 2Mj r . 

In section II we present the mathematical preliminaries of the TEGR, its Hamiltonian 
formulation and the expression of the energy for an arbitrary asymptoticaly flat spacetime. 
In section III we carry out the construction of triads for a three dimensional spacelike 
hypersurface of the Kerr type, obtain the general expression of the energy contained in 
a volume V of space and provide the exact expression of the latter in the limit of slow 
rotation. Comments and conclusions are presented on section IV. 

Notation: spacetime indices [i,is, ... and local Lorentz indices a, b, ... run from to 3. 
In the 3+1 decomposition latin indices from the middle of the alphabet indicate space 
indices according to \i = 0, i, a = (0),(i). The tetrad field e% and the spin con- 
nection uj^ ab yield the usual definitions of the torsion and curvature tensors: R a b^u = 
d^uj u a b + cOfj, a c ^v c b — T a ^ v = d^e a v + cu^ a be b v — .... The flat spacetime metric is 
fixed by r\ {m) = -1. 
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II. The TEGR in Hamiltonian form 



In the TEGR the tetrad field e a M and the spin connection are completely inde- 
pendent field variables. The latter is enforced to satisfy the condition of zero curvature. 
The Lagrangian density in empty spacetime is given by 

L{e,u,\) = -ke(-T abc T abc + ^T abc T bac - T a T a ) + e\ ab ^R abfll/ {u) . (1) 

where k = j^q, G is the gravitational constant; e = det(e a M ), \ ab ^ u are Lagrange 
multipliers and T a is the trace of the torsion tensor defined by T a = T b ba . 

The equivalence of the TEGR with Einstein's general relativity is based on the identity 

eR(e,u) = eR(e) + e(^T abc T abc + T abc T acb - T a T a ) - 2d^eT») , (2) 

which is obtained by just substituting the arbitrary spin connection u^ ab = °uj^ ab {e) + 
K^ ab in the scalar curvature tensor R(e,u>) in the left hand side; u^ a b(e) is the Levi- 
Civita connection and K^ ab = \e a x e b u (T\ flu + T u \^ — T^\) is the contorsion tensor. The 
vanishing of R a b/JiU (u), which is one of the field equations derived from (1), implies the 
equivalence of the scalar curvature R(e), constructed out of e a M only, and the quadratic 
combination of the torsion tensor. It also ensures that the field equation arising from the 
variation of L with respect to e a M is strictly equivalent to Einstein's equations in tetrad 
form. Let = denote the field equations satisfied by e a/i . It can be shown by explicit 
calculations that 

^7 = \{ R M ~ l^R(e)} . (3) 
(we refer the reader to ref.@ for additional details). 
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It is important to notice that for asymptotically flat spacetimes the total divergence 
in (2) does not contribute to the action integral. This term is a scalar density that falls 
off as ^3 when r —>■ oo. In this limit we should consider variations in g^ u or in e afl that 
preserve the asymptotic structure of the flat spacetime metric; the allowed coordinate 
transformations must be of the Poincare type. The variation of d^{eT^) at infinity un- 
der such variations of e aAt vanishes. Moreover all surface integrals arising from partial 
integration in the variation of the action integral vanish as well. Therefore the action 
does not require additional surface terms, as it is invariant under transformations that 
preserve the asymptotic structure of the field quantities. This property fixes the action 
integral, together with the requirement that the variation of the latter must yield Ein- 
stein's equations (the Hilbert-Einstein Lagrangian requires the addition of a surface term 
for the variation of the action to be well defined; a clear discussion of this point is given 
in ref. [|12| ). In what follows we will be interested in asymptoticaly flat spacetimes. 

The Hamiltonian formulation of the TEGR can be successfully implemented if we 
fix the gauge u 0a b = from the outset, since in this case the constraints (to be shown 
below) constitute a first class setO. The condition u>o a b = is achieved by breaking 
the local Lorentz symmetry of (1). We still make use of the residual time independent 
gauge symmetry to fix the usual time gauge condition e^)° = e( )i = 0. Because of 
^>oab = 0, H does not depend on p kab y the momentum canonically conjugated to u>kab- 
Therefore arbitrary variations of L = pq — H with respect to P kab yields iotab = 0. Thus 
in view of u 0a b = 0, u^ah drops out from our considerations. The above gauge fixing can 
be understood as the fixation of a global reference frame. 

Under the above gauge fixing the canonical action integral obtained from (1) becomes|| 




(4) 
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H = NC + N*C t + S mn n m " + J-d k (NeT k ) + d k {W k Nj) . (5) 

07TG 

iV and N l are the lapse and shift functions, and S mri = — E nm are Lagrange multipliers. 
The constraints are defined by 



C = dj(2keT j ) - keZ kij T kij - ^(n%* - ^IT 2 ) , (6) 



C k = -e^d^ - n^% )ik , (7) 

with e = det(eQ) k ) and T l = g lk e^ l T^y k . We remark that (4) and (5) are invariant under 
global SO (3) and general coordinate transformations. 

We assume the asymptotic behaviour e^) k ~ r}j k + \hj k {^) for r — > oo. In view of the 
relation 

^ | tfxdjieT*) = J s dS k (d t h tk -d k h u ) = E ADM (8) 

where the surface integral is evaluated for r — > oo, the integral form of the Hamiltonian 
constraint C = may be rewritten as 

J d 3 x{keX kij T kij + ^(n%* - ^n 2 ) | = E ADM . (9) 

The integration is over the whole three dimensional space. Given that dj{eT^) is a scalar 
density, from (7) and (8) we define the gravitational energy density enclosed by a volume 
V of the space as 



E > = ic Iv d3x8 ' {eT,) ■ (10) 
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It must be noted that E g depends only on the triads e^j restricted to a three-dimensional 
spacelike hypersurface; the inverse quantities e^ k ' 1 can be written in terms of e^. From 
the identity (3) we observe that the dynamics of the triads does not depend on u>^ a b- 
Therefore E g given above does not depend on the fixation of any gauge for ui^ a b- 



III. Energy of the Kerr geometry 



The Kerr solution|TB[ describes the field of a rotating black hole. In terms of Boyer 
and Lindquist coordinates JT^] (t, r, 8, <p) it is described by the metric 



ds 2 = -^-[dt-asin 2 9d(j)} 2 + ^-[(r 2 + a 2 )dcf) - adtf + ^dr 2 + p 2 d9 2 , (11) 
p 2 p 2 A 



A = r — 2mr + a , 

2 2,2 2n 

p = r + a cos 9 ; 

a is the specifc angular momentum defined by a = — . The components of the metric 
restricted to the three dimensional spacelike hypersurface are given by g n = c/ 2 2 = p 2 
and (?33 = ^-sin 2 9, where S is defined by 

T 2 = (r 2 + a 2 ) 2 - Aa 2 sm 2 9. 
We define the triads e^j as 



( -^=sin9 cose, 
-4=sin9 sine 

VA 1 

-ttCos9 



pcos9 cose, 
pcos9 sine, 
—psin9 
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——sin9 sin<j)\ 

— sin9 cosS 
p Y 
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(k) is the line index and i is the column index. The one form is defined by 



e <*) 



e (k) r dr + e (k) e d6 + e (fc) # 



from what follows 




= p ldr 2 + pW + Lm 2 ^ 2 



■4=sm#. Therefore the triads given by (12) describe the 



We also obtain e = det(e^)i) = 



components of the Kerr solution restricted to the three dimensional spacelike hypersurface. 

One readily notices that there is another set of triads that yields the Kerr solution, 
namely, the set which is diagonal and whose entries are given by the square roots of gu. 
This set is not appropriate for our purposes, and the reason can be understood even in the 
simple clase of flat spacetime. In the limit when both a and m go to zero (12) describes 
flat space: the curvature tensor and the torsion tensor vanish in this case. However, for 
the diagonal set of triads (again requiring a — > and m — > 0), 



some components of the torsion tensor do not vanish, T( 2 )i2 = 1, ^(3)13 = sin9, and E g 
calculated out of the diagonal set above diverges when integrated over the whole space. 
Therefore the use of (12) is mandatory in the present context. 

The components of the torsion tensor can be calculated in a straightforward way from 
(12). Only T( 3 ) 13 and T^ 2 3 are vanishing. The others are given by: 



e M = dr , e ie) =rd6 , e w 



r sinO dcp 



cosOcoscj) ( — h 



2 

r or 



P 



) 



(1)12 



p pVa 
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T (1)13 = sinOsinM — \- [2r(r 2 + a 2 ) - a 2 sin 2 6 (r - m)] + ^ + -^=} 

pE p^ VA 



£ 2 - 2^ A E 



7)i)23 = cosOsin(f){p V a sin 9 (— -)} 

P P 



2 

7)2112 = cos9sind>( — I -=sin 2 9 -=) 



T (2U3 = -sinOcosc/) { [2r(r 2 + a 2 ) - a 2 sin 2 (r - m)] + — + -^=} 

pS p 4 VA 



7)2)23 = -cos9cos(j) {p - — + a 2 sin 2 9 (-^ - -^)} 

p pE p d 
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T (3)12 = ^[-^ + -| + -^= C0 ^] 
In order to evaluate (9) we need to obtain T\ After a long calculation we arrive at 



Tl = VA + VA _ A [2r(r 2 + a2) _ a 2 sm 2^ (r _ m)] 
p z E p z E z 



m , . „ „ a 2 1 cos9 E 2 . 2 A E 

T = sin9 cos9 — + — — — [p + a 2 sin 2 9(— - — , 

p 4 pE sin 9 p ph p 4 



T A = . 

The gravitational energy density inside a volume V of a three dimensional spacelike 
hypersurface of the Kerr solution can now be easily calculated. It is given by 



E y^A 

sin9[p H I 2r(r 2 + a 2 ) — a 2 sin 2 9(r — m) 

p pE 
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d_ 

+ d6 



£a 2 9 „ „ cosO ( £ 9 9 „,A £ 
sm 2 9cos9 + 1 - - 



, p ha sm 9(— — - 

Ap 3 p pS p 3 



(13) 



Next we specialize E g to the case when the volume V is contained within a surface with 
constant radius r = r D assuming r a > r+, where r + 



ni 



a 2 is the outer horizon 



of the black hole. The integrations in and r are trivial. Also, because we integrate 9 
between and tt, the second line of the expression above vanishes. We then obtain 



1 r f E 

E„ = - d9sin6<p-\ 

4 Jo I p 



2r(r 2 + a 2 ) — a 2 sin 2 9(r — m) 



(14) 



p pS 

We have not managed to evaluate exactly the integral above. However, in the limit 
of slow rotation, namely, when — << 1 all integrals have a simple structure and we can 
obtain the approximate expression of E g . It reads 



E„ 



rj 1 



\ 



2m a 2 \ <r 
1 + + 



' 2m / 2m\ 



2m a 2 
1 + - 



(15) 



r K J 6r Q 

This is exactly the expression found by Martinez || for the energy inside the surface of 
constant radius r Q in a spacelike hypersurface of a Kerr black hole, in the limit of small 
specific angular momentum. As in ref.||, we have not expanded the square root which 
appears in (15) in powers in % . 

We remark that the expansion of p + — in the integrand of (14) yields — eo, whereas 
the remaining term corresponds exactly to e, expressions (3.17) and (3.1) respectively of 
0. It does not seem to be possible, however, to split di(eT l ) into two terms such that 
their integrals arise in the form e — 6q. 

As a very interesting application of (14), let us calculate the energy contained within 
the outer horizon, i.e., we will calculate (14) when the surface of constant radius is defined 
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by r = r + . This surface is characterized by A = 0. The integral can be calculated exactly 
for any value of a. The latter is parametrized in terms of the black hole mass m according 
to 



a = km , < k < 1 
After a number of integrations we arrive at 



Eg — m 



J2p~ 6p-k 2 , {J2p~+k 

+ T, l- n 



4 Ak V p 



(16) 



where p is defined by 



V 



1 + VT^A^ . 



This is the ammount of energy expected not to escape from the black hole by any process 
in which the black hole interacts with external particles. It is then important to compare 
(16) with E = 2M ir . 

We recall that a rotating black hole can have its mass decreased by means of Penrose's 



process of extraction of energy |[T5[|. The idea is the following. We consider a particle that 
is emitted towards the black hole and penetrates into the ergosphere. Suppose we arrange 
the particle to break up into two fragments, in such a way that one of the fragments has 
total negative energy. This is possible in principle, since the energy need not be positive 
in the ergosphere. By conservation of energy, the fragment with positive energy has an 
energy greater than that of the incident particle. Thus energy will be extracted from the 
black hole if the positive energy particle escapes to infinity, while the black hole absorbs 
the negative energy one. As a consequence, the mass of the black hole is decreased. We 
expect, however, that not the whole energy of the black hole can be extracted in this 
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manner. The existence of the horizon certainly prevents one from exhausting the total 
energy. ChristodoulouJTTJ has given an argument to determine how much energy can be 



extracted from the black hole by Penrose's process. He concluded that at the end of 
this process (when the ergosphere disappears and the black hole becomes static) the final 
(irreducible) mass of the black hole is given by 



1 



M ir = -\}r\ + a 2 . 

Martinez U has calculated the energy inside the horizon of the Kerr black hole in the 
limit of small a. He arrived at E = 2M ir [l + O(jjr)]. A similar result (in the same 

ir 

approximation) has been obtained by Zaslavskii||16|| in the analysis of a generic axially- 
symmetric spacetime. The question immediately arises as to whether this relationship 
holds for any value of a. This is in fact the conjecture made in ref.0. 

Since expressions (14) and (16) are valid for any value of a in the appropriate range, 
it is worth comparing (16) with 2M ir . In our parametrization we have 



E = 2M ir = m^2p . (17) 

The expression above certainly looks different from (16). However in the range < 
a < m expressions (16) and (17) as functions of k are strikingly indistinguishable, as we 
can see in Fig.l. In the latter we have plotted — against k. The upper curve represents 

(16) , the lower one (17). We see that for small values of the parameter k the two curves 
are essentially coincident. A tiny deviation occurs for values of k near 1. Inspite of this 
deviation, this is a remarkable result in favour of (14). 

Unfortunately we have not been able to explain such small deviation between (16) and 

(17) for values of k near 1, although we expect such explanation to be of fundamental 
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importance. It might be related to some physical property of the Kerr black hole which 
we do not understand yet. 

Comments 

The gravitational energy E g defined by (14) can be evaluated for an arbitrary value 
of a by means of numerical integration. This is the major advantage of our procedure 
as compared to that of Brown and York||. By means of the latter one cannot construct 
expressions like (13) and (14), which may be useful in the study of astrophysical prob- 
lems, since in a general situation Brown and York's procedure requires the embedding of 
an arbitrary two dimensional boundary surface of the Kerr space in the reference space 
E 3 , a construction which is not possible in general (the evaluation of Eq in || is only 
possible in the limit f- « 1). Therefore the present approach is more general than that 
of ref. 0. Finally we remark that we expect expression (10) to be useful in the study of 
the thermodynamics of self-gravitating systems, where the gravitational energy plays the 
role of the thermodynamical internal energy that is conjugate to the inverse temperature. 
We hope to come to this issue in the future. 
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Figure 1: ^ against k. The upper curve represents the energy expression (16), the lower 
one (17). 
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